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Complex scattering as canonical transformation: 
A semiclassical approach in Fock space 
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We show that a theory of complex scattering between 
many-body (Fock) states can be constructed such that 
its classical limit is a canonical transformation thus en¬ 
coding quantum interference in the semiclassical form 
of the associated unitary operator. Based on this idea, 
we study the different coherent effects expected under 
different choices of the many-body states and provide 
different representations of the associated transition 
probabilities. In this way, we derive exact relations and 
representations of the scattering process that can be 
used to attack timely problems related with Boson Sam¬ 
pling. 


1 Linear scattering of photons 

We consider a typical scattering scenario, where a highly 
coherent many photon state of light is injected through 
waveguides into a complex array of optical elements, 
such as e.g. in [1-5]. We further assume that decoher¬ 
ence and dephasing due to losses and /or coupling with 
uncontrolled degrees of freedom can be neglected. The 
simulation of such a scattering process between multi¬ 
particle photonic (or in general, bosonic) input and out¬ 
put states is a computationally hard problem because 
it involves, as shown bellow, the calculation of perma¬ 
nents of large matrices. The complexity of this problem 
is expected to render the problem of sampling the space 
of matrices with a distribution given by their perma¬ 
nents, the Boson Sampling (BS) problem, also hard. Thus, 
a quantum optical device that samples for us scatter¬ 
ing probabilities between many-body states constitutes 
a quantum computer that eventually beats any classical 
computer [6] in the BS task, an observation that has at¬ 
tracted enormous attention during the last years [7-11]. 

The physical operation of our scattering device con¬ 
sists of mapping the incoming many-photon states |in) 
into the output states |out). By injecting the same incom¬ 
ing state several times and counting the number of times 


we get |out) as output, we will eventually obtain the tran¬ 
sition probability 

P(|in> — |out)) := |A(|in) — |out))| 2 = |(out|in)| 2 , (1) 

and our goal is to study this quantity. 

As any other quantum state of the field, the |in), |out) 
states belong to the Hilbert (Fock) space J6 of the system, 
which consists of all possible linear combinations of Fock 
states [12] 


|n) := \ m,nz,...,nM) 


( 2 ) 


specifying the set of integer occupation numbers ri\ ./%/■ 

An occupation number n, specifies how many photons 
(bosons) occupy the zth single-particle state. The choice 
of these channels (or orbitals) is a matter of convenience, 
depending on the particular features of the system. In 
the scattering problem there are two preferred options 
to construct the Fock space, namely, by defining occu¬ 
pation numbers specifying how many photons occupy 
a given single-particle state with either incoming or out¬ 
going boundary conditions in the asymptotic region far 
away from the scatterer. The operators that create a par¬ 
ticle in the case of given incoming boundary conditions 
are denoted by /r , and their action on the vacuum state 
|0) produces Fock states in the incoming modes [12]: 


i„in\ . i „in ,,in in, r~r 

In ) := |«! ,n 2 = [ 


b I 


t]"i 


-| 0 >. 


(3) 


Any operator acting in can be written as a multilinear 
combination of the creation operators and their adjoints 
b, called annihilation operators. The operator algebra is 
thus uniquely fixed by the canonical commutation rela¬ 
tions 


[b it bj] = 0 and [b it b\] = 8 ij . 


(4) 
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Similarly, the operators d f create photons in the single¬ 
particle states defined by outgoing boundary conditions, 
and representing physically photons exiting the scatter¬ 
ing region along a given channel. A fundamental obser¬ 
vation is that the Fock space can be equally well con¬ 
structed out of the many-body states defined by specify¬ 
ing occupation numbers in the single-particle outgoing 
states: 


|n out ) := |nr 1 , n° ut , ..., n™ x ) = Y\ |0>. (5) 

«■ 

The relation between incoming and outgoing Fock states 
is fully determined by a single-particle property, namely, 
the transition amplitude of the single-particle process 

o u = <0,...,0,n° ut = l,0,...,0|0,...,0,nf = 1,0,...,0), (6) 

which defines the single-particle scattering matrix with 
entries a;j. By comparison with Eqns. (3) and (4), we 
find 

dj = cs jib j (7) 

i 

which then allows us to relate the expansion coefficients 
cjf and c™*, appearing in the "in” and "out" representa¬ 
tions of an arbitrary many-body state, 

|t//> = XCln in > = LOm 0llt >, (8) 

n m 

through the amplitude 

A F (n,m) := (m out |n in ). (9) 

So far we have focused on the transformation proper¬ 
ties between Fock states, but the same questions can be 
addressed for other type of many-body states. Consider 
for example the common eigenstates of the incoming cre¬ 
ation operators [12] in the incoming basis, 

h,.|0 in > = </>f|0 in >, (10) 

so-called coherent states, which are labeled by a continu¬ 
ous set of complex numbers t />;. Although coherent states 
are not eigenstates of a commuting set of hermitian op¬ 
erators, they can be experimentally prepared [13], and in 
some sense they are the most classical states of the elec¬ 
tromagnetic field. Again, it can also be shown that both, 
in- and out-going coherent states, are an (over)complete 
basis of the Fock space, and the amplitudes 

A c i X ,(p)-=((f> 0Ut \X in ) (ID 


are the matrix elements of a many-body unitary transfor¬ 
mation performing the change of representation from in¬ 
coming to outgoing coherent states. 

The third basis set that we are going to discuss is 
defined by the common eigenstates |q m,out ) of the so- 
called quadrature operators [14] that correspond to the 
quantum operator associated with the observable elec¬ 
tric field [15] 

Q? := + b\ , <?° ut := d t + d ]. (12) 


It is easy to show that quadrature eigenstates 


-in,out I in,out 


.in.out, in.out 


|qin,oui> = ^,uuL| ql 


(13) 


are labeled by a continuous set of real variables and that 
they are normalized (to the Dirac delta), complete and 
orthogonal. We can again define the corresponding trans¬ 
mission amplitude 

7l Q (q, Q) := <Q° ut |q in >- (14) 


The construction of the transformations between the 
different basis sets is cumbersome but straightforward 
and we refer the reader to the references [12, 14] for 
further details. We just note that for a given choice 
of single-particle orbitals, all operators (number, cre¬ 
ation/ destruction and quadratures) commute with each 
other if they correspond to different single-particle states 
(or index i ). Therefore the results for a given mode [12,14] 
are sufficient, 


(q\n) = Jl_ -=H 

\j2 n n\\/2n 



(q\cp) = 


(2tt) 


1/4 


I0r 

2 


(f-</>) + 


number to quadrature, 
coherent to quadrature, 


(15) 

1 2/ 

{n\(p) =— =(f) n e^^ /2 , coherent to number, 

Vn\ 

where H n {q) is the n - 111 Flermite polynomial. 


2 The Boson Sampling problem 

2.1 Outline of the problem 

With the formalism presented in the last section, we re¬ 
turn to our original problem, namely the explicit calcula¬ 
tion of scattering amplitudes between Fock states. Using 
Eq. (9) and the definitions in Eqns. (3) and (5) we get the 
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exact expression 


A f (n, m) := <0| 



| 0 >. 


(16) 


Our goal is to obtain an expression for A v , and eventu¬ 
ally for the transition probabilities |A F | 2 , in terms of the 
single-particle scattering matrix a. At first glance, due to 
the absence of interactions this seems to be an easy task 
since in the scattering process the total amplitude fac¬ 
torizes in terms of the amplitudes of individual, single¬ 
particle processes. However, while this is the case for 
systems of non-interacting distinguishable particles, for 
the case of interest here quantum effects due to indis- 
tinguishability render the calculation of scattering ampli¬ 
tudes a hard problem [6,11], as it is apparent when try¬ 
ing to calculate the amplitudes by substitution of Eq. (7) 
into Eq. (16) and further using the commutator in Eq. (4) 
that the complexity of the result is combinatorial in ori¬ 
gin. Since the explicit calculation has been reported else¬ 
where, we just present the final result 


A f (n, m) = Perm M(<r). 


(17) 


and refer the reader to [6,11] for further details. Follow¬ 
ing [6,11], the transition amplitude obtained by this pro¬ 
cedure is given by summing up products of entries of a, 
where each term is actually a permutation of the multi 
dimensional indices labeling output channels. It can be 
therefore written in terms of a new matrix M(tr) (ob¬ 
tained by repeating the y'-th row of a n j times and the 
j -th column rri j times 1 ). The key observation is that we 
indeed sum all the terms obtained by permuting the sec¬ 
ond index of this enlarged matrix, resulting in an object 
known as permanent. 

In this way, the physical scattering of photons pro¬ 
vides a physical device that calculates permanents of 
large matrices. One only needs to repeatedly measure the 
output state and the accuracy with which the device cal¬ 
culates (or simulates) the precise value of the associated 
permanent can be made arbitrarily large by repeating the 
measurement as many times as needed. In a further step, 
by randomly changing tr the device can be used to sam¬ 
ple the space of matrices with a weight given by their per¬ 


1 The precise and primary definition of M is Mj ,t = (TdjWd k 0a)< 

with d(n) being the Af-dimensional vector defined by dj (n) = 
min{fce {1,...,M} : »jt-i < j < «*}. By permuting columns 
and rows of M one arrives at the definition above. 


manent. It is this task, the BS problem, where under cer¬ 
tain conditions it is expected that the quantum device 
beats any classical computer [6]. 


2.2 Many-body scattering as canonical 
transformation 

Small scattering devices calculating permanents can be 
actually realized, with several examples now available 
[10], while preparation techniques that allow for coher¬ 
ent creation of correlated photons beyond N - 30, where 
the sampling using the quantum device will beat any 
classical computer, is presently matter of intense re¬ 
search [6]. Realistic scenarios, however, seem to reach se¬ 
vere complications already for N - 12. Thus, it seems im¬ 
portant to explore whether the fundamental aspects of 
complex many-body scattering allow for other types of 
implementation, different from the photonic ones. Here 
we try to approach this question from a more abstract 
perspective. 

We only demand the matrix a to be a unitary single¬ 
particle scattering matrix, 

[O' -1 ];,; = [a-]*j. (18) 

This means that the calculation of permanents of large 
matrices is realized by a device with outcomes given by 
the amplitudes 

A F (n,m) := <m'|n>. (19) 

We call |m) the state specified by occupations m in the 
unprimed basis and |m') the state specified by m in the 
primed basis. The later is constructed out of the opera¬ 
tors 

h 'i = L u jA and ( h 'if = I u h fof (20) 

i i 

for any unitary matrix u. Note that for the choice u = a 
we recover the scattering version with b'j = dj. 

Finally, a straight forward calculation shows that, 

[b' i} b'j] = 0 and [£',(£')}] =S i} (21) 

follows from Eq. (4). The transformations (20) are linear 
and canonical, the latter because they do not change the 
algebraic relations between the basic operators. 

We conclude that the BS problem can be realized by 
any nontrivial device where transition amplitudes are 
measured between Fock states defined by two different 
sets of creation operators, defined with respect two dif¬ 
ferent single-particle basis sets. The physical implemen¬ 
tation of BS requires then, first of all, a measurement 
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protocol that provides the many-body transformation be¬ 
tween Fock states after a linear canonical transforma¬ 
tion. 

Note that for the BS problem, the essential property of 
the scattering device is that it provides permanents, and 
the possibility of connecting permanents with transition 
amplitudes is entirely due to the linearity of the single¬ 
particle transformation. Any physical system where the 
mapping b —■ b' is nonlinear, as happens for general 
(non-linear) unitary transformations 

b-+b'=U f bU, U = e iG(bb+) , (22) 

with a hermitian but non-quadratic generator G, defines 
a quantum device that still calculates transition ampli¬ 
tudes but not permanents. 

Using this broader view, it is in principle possible to 
implement other processes where the output of a mea¬ 
surement is given by permanents and therefore can be 
used as a basis for the physical implementation of BS. 
The quantum optical scenario involving scattering of 
photon states has some very attractive features, in par¬ 
ticular that the many-body output states can be indeed 
measured at the single-photon level by photocounting, 
while its main drawback is the difficulty to prepare of 
photonic Fock states with large total number of photons 
(the state of the art is N = 6). Since quantum states of 
indistinguishable bosonic atoms with macroscopic oc¬ 
cupations can be prepared by cooling techniques, cold- 
atoms alternative offer an interesting possibility for BS. 
The drawback here is the difficulty of performing tomog¬ 
raphy of many-body cold atom systems at the single¬ 
particle level, but advances in this direction are under 
way [16]. 

Assuming for the moment that the measurement 
of many-body states in cold atom systems reaches the 
regime of single-atom precision, we sketch a possible BS 
scenario for such systems. Consider a system of ultracold 
atoms in an optical lattice, where the hopping amplitude 
between adjacent sites is / and the strength of the in¬ 
terparticle interaction is V. Assume now that initially (at 
time t~) we have V » /, and the interaction energy is so 
large that hopping gets completely suppressed [17], the 
so-called Mott phase. In good approximation the ground 
state of the many-body system is a Fock state where the 
occupations refer to the number of atoms in each site, 
namely, a Fock state constructed out of single-particle 
states defined by localized (Wannier) orbitals [18], The 
"quench" scenario is defined by an abrupt change of pa¬ 
rameters (possible by tuning the atoms through a Fesh- 
bach resonance [19,20]) at time t + , such that we have 
now / » V. We are interested in the transition ampli¬ 
tudes between the initial state and the eigenstates of the 


quenched Hamiltonian, where the later are again Fock 
states but built from delocalized (momentum) single¬ 
particle orbitals. The calculation of these transition am¬ 
plitudes is strictly given by A F (n,m), with the specific 
choice for u as the matrix that linearly relates the Wan¬ 
nier and the momentum orbitals. Thus, such a device 
provides permanents as its output. Furthemore, if the on¬ 
site energies in the Mott phase are chosen to be random, 
the matrix u is itself random, and BS can be fully imple¬ 
mented. 


3 Equivalent representations 

We return now to the question of how a general single¬ 
particle linear canonical transformation is reflected in 
the transformation of the different (Fock, quadrature and 
coherent) many-body states and how the hardness of cal¬ 
culating permanents gets reflected in the different repre¬ 
sentations. 


3.1 Coherent states 

The simplest transformation between many-body states 
after a single-particle canonical transformation is the 
one for the coherent states, and hence we start with this 
case. Any coherent state can be constructed out of the 
vacuum state |0) by the application of the displacement 
operator [12] 

D(^,^*) = e^*' 6+< ^' 61 (23) 

as 

I0> = D(<p,(f)*)\0), (24) 

and similarly for the primed states | y/') 

\y/') = e v '* b ' +vb ' + |0>. (25) 

From this, and the defining relation between primed and 
unprimed canonical operators in Eq. (20), we get 

\y/') = \uy/). (26) 

Using again well known properties of the coherent states 
[12], the transition amplitude is given by 

A C ((P, i/r) = e - ***’■*-W-V+V* (27) 

This result implies in turn for the corresponding transi¬ 
tion probability 

p c ((p,y/) := \A c {(f>,y/)\ 2 = (28) 
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admitting a straightforward interpretation, very much 
consistent with the idea that quantum coherent states 
are the most classical states of light: At the classical level, 
the canonical transformation simply consists of a linear 
transformation between the held amplitudes given by 
(p —► u • (p. The classical probability to obtain the state y/ 
after a canonical transformation of the state (p is imple¬ 
mented is nonzero only if y/ = u • (p. In the quantum case, 
this sharp peak is smoothed into a Gaussian. 

In terms of the scattering scenario, the transition 
probability between coherent states also agrees with in¬ 
tuition: The probability is strongly peaked around the 
output state labeled by the classical held amplitude re¬ 
sulting from scattering of the classical input held. 


3.2 Quadrature states 

In the same spirit as in the case of coherent states, 
the transformation rule for quadrature states can be de¬ 
duced by the corresponding transformation for the defin¬ 
ing canonical operators. In the coherent state case, the 
canonical pair is b, if , and therefore for the quadrature 
case we must hnd the set of canonical conjugate partners 
of the q’s. The obvious choice that turns out to do the job, 
is to dehne [14] 


<7,•:=£,. +fit , 

As the rp quadratures, the p-quadratures have a com¬ 
plete, orthogonal and Dirac-normalized common set of 
eigenstates, 


Pi lp> = Pi lp>. 


(29) 


The analogy with the usual position and momentum op¬ 
erators in particle (hrst quantized) quantum mechanics 
is evident after using their definition to obtain [14] 


<qlp> = 


eiw 

(47l) M/2 


(30) 


However, it must be stressed that quadrature states do 
not represent any single-particle property at all. In fact, 
it can be shown that they do not represent states with a 
well defined total number of particles, thus making their 
interpretation as any sort of localization property in real 
space impossible. 

Our goal is again to inter-relate the two quadrature 
states IQ') and |q), defined by 

qlq> = qlq>, (31) 

q'IQ'> = QIQ>, 


using as input the canonical transformation given by 

q+ ip —■ q' + ip' = u(q+ ip). (32) 

This canonical transformation can be solved for q' sim¬ 
ply by taking its hermitian part on both sides using the 
decomposition 

u = u r + in' (33) 


into real and imaginary parts. The eigenvalue equation 
defining |Q'> is then found to be 


i d 1 , _ . 

/U 1. ( u l. q _Q) 


dq 2 


<Q'lq> = 0. 


(34) 


This can be solved using a Gaussian ansatz to get 


A Q (q,Q):=< Q|q> 



with similar expressions for the p-quadrature states. 

Using Eq. (35), we obtain an interesting result for the 
transition probability between quadratures, 

P Q {q,Q):=\A^{q,Q)\ 2 =- --- — (36) 

det47Tu(u ! ) 


It is fully independent of the initial and final states. In 
the scattering scenario this means that the probability to 
obtain a given configuration after measuring the electric 
field in the output channels is the same for any input and 
output configuration. As it is clearly seen in Eq. (35), how¬ 
ever, the amplitudes themselves are very structured func¬ 
tions of the input and output quadrature states and it is 
only the associated probabilities that display a flat pro¬ 
file. 


4 Exact representations 

Armed with the results of the last section we can now 
construct different exact expressions for the transition 
amplitudes between Fock states, A F (n,m), that, supple¬ 
mented with an ensemble of random u matrices, provide 
different representations of BS. Since the transition am¬ 
plitudes, Eqs. (27,35), in both quadrature and coherent 
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representation are not difficult to evaluate, the complex¬ 
ity of calculating permanents must stem from the trans¬ 
formations between the different basis sets. In the follow¬ 
ing we will make this connection explicit. 

Using the transformation rules, Eq. (15), between co¬ 
herent, quadrature and Fock states we obtain [12,14], 


A F (n,m):= <m'|n> = ^7 J dy/d(p{m'\y/')(y/'\<t>)((p\n) 

r ^w mi [(b*) ni e-^ 2 - 3 ^« |2 „ 

= J diffdtp FI—- ! _ 9 7 = = , - A c {<p,iff), 


n 2 y/milrii\ 


and 


(37) 


(n,m) := <m'|n> = J dQdq{m'\Q')(Q'\q){q\n) 


r H -»t [zk H "t § e “^“^ n 

= j dQdq fl- 1 J --f " - A Q (q,Q). 


v 2 n ‘ +m ‘ +1 7tn i \m i \ 


(38) 


Equations (37) and(38) are two equivalent representa¬ 
tions of the scattering amplitudes and provide a basis for 
realizing BS when supplemented with a physical ensem¬ 
ble of unitary matrices u. The first expression in terms 
of the coherent state transition amplitude A c is conve¬ 
nient for exact calculations, while the second equation in 
terms of A (> will be important when we connect BS with 
a three-step canonical transformation in order to under¬ 
stand its asymptotics for large N. 


5 A generating function for transition 
amplitudes 


sion, 


A (n, m) = 


Vmtlriil dx™‘ dy"‘ 


n 

i 

n s/mi'.rii'. 

(V t-4* 2 ) 


grrii grii 




e xuy 


x=y=0 


f i—r dxjdyi 
f m, + l n, + l 

J 1 X i Vi 


(39) 


e x " y 


which is one of our main results. It is a generating func¬ 
tion providing the transition amplitudes as high-order 
derivatives of an multivariate exponential function and 
generalizes [21]. Here it is clear that in any representa¬ 
tion the complexity of many-body scattering comes from 
the combinatorics involved in taking high-order deriva¬ 
tives of large products of exponentials. The second ex¬ 
pression, obtained by using the Cauchy integral formula 
(the closed integration contours enclose the origin), fur¬ 
ther transforms the problem in a way suitable for asymp¬ 
totic analysis. 

The generating function approach provides a way to 
eventually address some open questions, in particular 
the calculation of high order moments of the distribu¬ 
tion of transition amplitudes (or transition probabilities) 
over the ensemble of single-particle canonical transfor¬ 
mations [6], The particular advantage of this representa¬ 
tion is that the average over the unitary group of matri¬ 
ces u representing the single-particle canonical transfor¬ 
mation can be performed exactly. In section (6) we show 
how to follow this program in the simpler case of Ginibre 
(complex) matrices, and provide for the first time exact, 
explicit expressions for the third moments of the distri¬ 
bution of squared permanents. 

So far, all equivalent versions of the transition ampli¬ 
tudes have been obtained by exact, identical transforma¬ 
tions. In the rest of this section we will focus on the partic¬ 
ular regime of high densities, i.e, when AT := Z/ n; » M, 
where we can safely assume that the majority of configu¬ 
rations satisfy 


It is instructive to show how one finds yet another version 
of the transition amplitudes using the equivalence of the 
two representations in Eqns. (37) and (38). To this end we 
use the identities 


n d ’ 1 

cb = -e 

r dk n 


k(p 


k=0 


, H n [q) = e c < 2 -^e- ( ^ fc)2 
1 dk n 


k =0 


which allow us to perform exactly the integrals over the 
intermediate variables y/,(p in the coherent state repre¬ 
sentation and Q, q in the quadrature case. After some 
calculations we get the exact, surprisingly simple expres- 


rij » 1, nq » 1, (40) 

and powerful methods of asymptotic analysis can be 
safely applied. However, although BS involves the regime 
of large N and M, it is expected to be a hard problem 
only in a specific asymptotic limit given by M » N 2 [6], 
and the high-density limit cannot be used to make state¬ 
ments about it. The study of the behavior of the scatter¬ 
ing amplitudes in the appropriate dilute limit of interest 
for BS is currently under progress. 

If the conditions in Eq. (40) hold, we can then eval¬ 
uate the contour integrals in Eq. (39) by the method of 
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steepest descent applied to 


A b (n, m) = 


n 

i 


Vmi'.rii'. 


-4n 2 


f 


(41) 


\\dxidyi 

i 


3 -Zi(»i + l)logn-Zi(mi + l)logyi+X'Ky 


thus making contact with the theory developed in [22]. 
Here we are not interested in the technical details of the 
full calculation of the large-AT asymptotics, but instead in 
the physical interpretation of the saddle point conditions 


d 

dxi 


- X n, log Xi - X nit logy; + x u y 
i i 


= o, 


(42) 


d 

dyi 


~Y J n i log*/ - X m i logy; + X- u ■ y 

i i 


= 0 


(43) 


selecting the optimal values of the, so far, purely formal 
complex variables x,y. Under the variable transforma¬ 
tion 


Xi = y/nje l6 ‘,yi = s/rrqe ’ 11 


(44) 


The semiclassical approximation for the amplitudes that 
define the unitary operators representing the first and 
last canonical transformations of the chain in Eq. (47), 

7l qn (n,q) = (q|n> , H qn (N,Q) = <Q|N>, (50) 


is given by [23,24] 


N 


A qn (n,q)-r 


1 d 2 f{rii, qj) j f ( njilli ) 


=i V driidqi 


A qn (N,Q)~ Y 


~ 1 d 2 F(Ni,Qj) 


(51) 


,iF(N iy Qi) 


f = \ V 2ni dNidQi 


in terms of generating functions f{n,q),F[N, Q) = f(N, Q) 
satisfying 


e =rJ^’ Q =ik fiN ’ Q) - 


(52) 


Finding these generating functions is a standard problem 
with explicit solution 


f(n,q ) = -j^/4 n- q 2 - n arccos 


\2 y/nj 


(53) 


the resulting set of 2 M complex equations can be re¬ 
duced to find the M real angles xi satisfying the condi¬ 
tions 

X u i iu ir y/mim,,e l[xl ~ Xl ' ) = rii. (45) 

In other words, the asymptotic limit of many-body tran¬ 
sition amplitudes for large densities is dominated by con¬ 
figurations x, y satisfying 


Interestingly, and contrary to /\ qn and 7lQ N , due to the lin¬ 
earity of the transformation (q,p) — (Q,P) the interme¬ 
diate step (q,p) — (Q,P) (responsible for the change in 
single-particle representation) is not only approximated 
but it is in fact exactly given by the semiclassical expres¬ 
sion. The result is then identical to H^(q,Q) in Eq. (35). 

We can now construct the semiclassical approxima¬ 
tion for the full transformation (n, 6 ) —* (N, 0) by operator 
multiplication of the three intermediate transformations, 


y= ux with |y/| 2 = m; and |x,| 2 = (46) 

This shows that in the limit of large densities, the calcula¬ 
tion of transition amplitudes requires the solution of (46), 
namely the calculation of the phases of the classical in¬ 
put and output field amplitudes (linearly related through 
u) required to satisfy shooting (instead of initial-value) 
boundary conditions. This interpretation can be made 
even more explicit by considering the quadrature repre¬ 
sentation of the amplitudes. To this end, we consider the 
chain 


7l F (n,N) = JrfqrfQ(H QN (N,Q))*H Q (Q,q)^ qn (n,q), (54) 
to get 



(n, 0) —► (q,p) —► (Q,P) —► (N,0) 
defined by 

qt + ipt = y/rqe l9i , Q, + iP t = -/wie' 0 * 
and 

Q+ z'P = u(q+ ip). 


(47) 


(48) 


(49) 


n- 


expt i 


- X V 4 N J~Qj + N i arccos 


2ni\!Nj -Qjl4 


*n 

j 


exp 


| i ^ yj4rij - q 2 -nj arccos 


2 V W J 


IniJuj - q 2 14 
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(55) in order to write 


This is exactly the same result we obtain by considering 
the large n limit of the exact representation, Eq. (38), by 
using the asymptotics 


H n {q) = 


\ 


2n+l fliia-n+q 2 


1- 


r 


2n +1 


X cos 


n + - arcsin 
2 


_. -I— \ /2/t + 1 — cj~ 

V2^TT) 2 V 



(56) 


Note that the complexity of many-body scattering is re¬ 
flected in the coherent sums over quantum mechani¬ 
cal amplitudes explicitly appearing in Eq. (17), namely, 
quantum interference results in the highly irregular pat¬ 
tern one obtains for the transition probabilities as a func¬ 
tion of the incoming and output states [11]. Very much 
opposite to the semiclassical method presented here, 
quasiclassical approaches, based on adding probabilities 
instead of amplitudes, capture only the gross features of 
these patterns. To stress this point, it is important to un¬ 
derstand where quantum interference is hidden in our 
semiclassical approach. In terms of Eq. (17), by expand¬ 
ing the permanents of M(<r) as sums over products of 
single-particle scattering matrices, these coherent sums 
over products of single-particle paths can be made very 
explicit, as in [25], 

The semiclassical interpretation of many-body scat¬ 
tering (at least for the case of large occupations) allows 
us to understand the complexity of the problem, and the 
origin of massive quantum interference in terms of classi¬ 
cal canonical transformations. To this end, consider now 
the unique canonical transformation implementing the 
full change of canonical variables (n, 0) —<■ (N, 0) without 
the intermediate steps in terms of quadratures. Then the 
semiclassical theory of quantum canonical transforma¬ 
tions indicates that we must find the generating function 
ic(n,N) that, together with the definitions 


8=—w{ n,N) 

on 

, &=—w{ n,N), 
dN 

(57) 

V'tye 40 ' 

= L««7v TTJe"^, 

(58) 


i 


gives the explicit form of the transformation as 
N = N(n,0) , 0 = 0(n,0), (59) 


A F (n,N)oc 


det 


d 2 u;(n,N) 

dndN 


e i if(n,N) 


(60) 


However, in this case we encounter a new issue that was 
not present in the canonical transformations we have 
seen before: although the initial value problem of find¬ 
ing (N, 0) from (n, 0) admits a unique solution (given by 
the transformation equations), the boundary problem of 
finding (0,0) for given (n, N) admits a very large set of so¬ 
lutions. Each of these solutions represents a branch y of 
the multi-valued generating function w, and the correct 
form of the semiclassical approximation to the transition 
amplitude is then, 


1 

2 

e iu/ r tn,N)+iju r f ^ jgjj 

Here the index p r is a topological property of the particu¬ 
lar branch that can be computed from the classical trans¬ 
formation. As expected, this is also the solution of the 
calculation of the amplitudes using the generating func¬ 
tion (55), within the saddle point approximation. Hence 
the semiclassical origin of both, the complexity of many- 
body scattering and the massive quantum interference 
associated with it, is the highly non-linear form (and 
therefore the multi-valuedness) of the boundary prob¬ 
lem connecting occupations. 


A f ( n,N) = £ 

r 


det 


1 d 2 ic r (n,N) 


2 it dndN 


6 Distribution of permanents 


In this section we will calculate the first three moments of 
the distribution of permanents over the (complex) Gini- 
bre ensemble to show exemplary how the representa¬ 
tion (39) leads to a solvable combinatorial problem. The 
calculation is exact, in that it does not involve any asymp¬ 
totics. It would be of course important to perform a simi¬ 
lar calculation in the regime of interest for BS, and this is 
work under progress. 

Let cr 2 denote the variance of the independent real 
parts and imaginary parts of all matrix elements in A and 
let N be its dimension. We start with an exact representa¬ 
tion obtained from (39), in a slightly different form 


Perm A 


\ 

n 


\i =1 


d 2 ' 
dx i dy ii 


e xAy 


x=y=0 


(62) 


where x= (jtq,..., x,y) (and similarly for y) is a column vec¬ 
tor and r is transposition. This representation allows the 
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Gaussian average to be performed exactly. Define the ten¬ 
sor 


p (fc) = 



(63) 


such that 



= Tr 


2 n 

a Ip® 


(64) 


The average of |Perm A| 2 ” is evaluated by separating real 
and imaginary parts of the matrix A to get 


(|Perm A| 2 ”) = 


' n N q2 q2 

,E Pi dxf^dyf™ dxf k W- k) 


x e 


2a 2 T?j =1 Ll l=1 


(2k-l) ( 2 fc-l) (21) (21) 
c. y Xj y. 


x=y=0 


which is equivalent to 


(65) 


2 n N 

(|PermA| 2 ") = coefficient of ]~[ H x®y® in 

fc=ii=i 

nil(l + 2 

i,j k,l 


( 66 ) 


The evaluation of the coefficients in (66) is related to the 
following combinatorial problem. First of all we can re¬ 
move the factor 2a 2 in ( 66 ) and in return eventually mul¬ 
tiply the overall coefficient with (2<r 2 ) nN . For each value 
of i = 1,..., AZall the x) k \ k = 1,...,2 n, have to appear ex¬ 
actly once. They come in pairs x[ fc 5 x ^' with k' odd and 
Z' even. We start by counting the number of ways to com¬ 
bine different factors in f|" ;=1 (l + x[ 2 fc_ 1 ) x[ 2() ) to get each 
variable (for fixed i ) exactly once. This is equivalent to 
counting pairings between n (representing the even in¬ 
dexes) and n (representing the odd indexes), which itself 
is equivalent to counting permutations of n. We write 
(p(D p\ 2 ) pm) or abreviated (P(i) p( 2 ) ••• pm) to adress 
a specific permutation P e S n . Specific pairs shall be de¬ 
noted by the corresponding column (^) = ( p f /c j). For all x- 
variables one has to count N independent permutations 
of n. Writing those one below the other will be referred to 
as table. 

For the y-variables again N permutations of n have 
to be counted. Since they come in combination with the 
x-variables in ( 66 ) they are not independent from the x- 
pairings. Each tuple (i, k, l) representing a pair in the N 
x-pairings actually comes with a fourth entry as a four¬ 
tuple ( i,j,lc,l ). This means that the pairs (/c, /) building 
up the y-pairings have to be taken from the x-pairings. 


In other words the y-pairings have to be a rearrangement 
of the x-pairings, keeping the (fc, Z)-indexes of all pairs. 
We will refer to this as a vertical rearrangement or per¬ 
mutation, depending on the context. In the process of 
rearranging identical pairs have to be taken distinguish¬ 
able (e.g. vertically swapping two identical pairs in the 
y-table has to be counted additionally) since the set of 
touples {Ui, ji,k,l),[i 2 ,j 2 >k,l)} is different from the set 
{(it, J 2 , k, l), (hijhk, l)} (if z'i / z 2 .Ji / j' 2 ) although the 
(/c, Z)-indexes of the two pairs involved are the same. 

(i) n = 1: There is trivially only one permutation for 
each Z concerning the x-variables. The same holds for the 
y-variables but there are N\ ways to vertically rearrange 
all the (})-pairs. We get 

(|Perm A| 2 > = ( 2a 2 ) N N\. (67) 

(ii) n = 2: The two different permutations of n = 2 are 

Z J i = (] 1 ) and Z J 2 = Q 2 ), which are incompatible, mean¬ 
ing they do not share any pair. Let N\ (Mi) and de¬ 

note the multiplicities of Pi and P 2 in the x(y)-table. The 
incompatibility implies M\ = N\,M 2 = N 2 . The number 
of ways to distribute these permutations on N twice is 
( NgNg ) 2 an< ^ ^ le num t> er vertical permutations of pairs 
is (ATi!) 2 (AT 2 !) 2 . We get 

N 

(IPermAI 4 ) = (2 a 2 ) 2N £ 8 LNa , N {M) 2 

N lf N 2 =0 

= (2ct 2 ) 2 W AZ!(AZ+1)!. ( 68 ) 

(iii) n = 3: The 3! = 6 permutations of n = 3 are 

(Pi,.... Pq) = ((1 2 3), (2 1 3), (1 3 2), (3 2 1), (2 3 1), (3 1 2)). 
Again we let N a and M a (a = 1,_,9) denote the multi¬ 

plicities of the permutations P a in the x- and y-table re¬ 
spectively. We define the 9 pair-counters p a 

Pi = Ni + N 3 , p 2 = N 2 + N 5 , p 3 = IV 4 + AT 6 , 

P 4 = N 2 + Ng, p 3 = ZVi + A/ 4 , p 3 = ZV 3 + AT 5 , (69) 

Pi — ZV 4 + ZV 5 , pg = N 3 + ZV@, pg = ZVi + n 2 

forthepairs(J),(i),(i),( 2 ),( 2 ),(§),(3),(3),(l) (inthat 
order). Taking into account (a) the multinomials for the 
distributions of the permutations among the N rows for 
both x and y, (b) the restriction to y-tables that are verti¬ 
cal rearrangements of the x-table and (c) the vertical per¬ 
mutation of identical pairs for y yields 

<|PermA| 6 ) = 


6 

N 

6 

( N 

(2n 2 ) 3jv n 

L 

5 I. Na.N n 

L 

a= 1 

(Na=0) 

<3=1 

\Ma= 0 ) 


x n <V(N),Pa(M) n 7T7 ,, ]lPa(M)!. 

a= 1 ltnUV a .ivj n .; a 
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(70) 

The 9 x 6 -matrix ^^| N) j has rank 5 so there are 5 in¬ 
dependent restrictions from ria<5p a (N),p a (M)- Also the re¬ 
striction Y.„ M a = N is contained when Y. a N a = N ap¬ 
plies. Thus (70) can also be expressed containing only 
6 sums. In the following form the number of sums is re¬ 
duced to 7, keeping one restriction and M\ independent. 


(|Perm A| 6 ) = (2 ct 2 ) 3JV (A/!) 2 f 


6 

1 

<3=1 


N 

L 

\N a =Oj 


N 

S I.N a ,N Y 
Mi=0 


n«Pa(N)! 


where the M a (a > 1) are given by 


(71) 


AT 2 — N\ + IV 2 ~ Mi , M 3 — N\ + A /3 — M \, 

M 4 = Ni + N 4 - Mi , M 5 = N 5 - Ni + Mi , (72) 


Mq — A/g — Ni + Mi 


and ( := 0 for k e N\{0}. Applying (71) we evalu¬ 

ate the scaled third moment (|Perm A| 6 )/(2cr 2 ) 3iV /(A/!) 3 
for the lowest N to 6, 18, 79, 140, ^2, 

cyq 276442 228754 3697434 48374363 12084328 55026632 

5536562^ 2903^ 7 Q 5 139 3^^39326 17 ^545903l8° 3 1562468^26 

33^f^§263 95^^612809| 2 %^710705§i 1 %^1996636^ r^ 84 Ar 
734825 ’ 1527701175 ’ 8793675 ’ 6186609 lU1 

1,..., 23 respectively and use this to estimate an assymp- 


totically exponential (as opposed to factorial) scaling of 
this quantity proportional to e AN N v (1 + G’ (^)) w ith A ~ 
0.3. 


7 Conclusions 

We have shown that the usual many body scattering 
scenario realizing the Boson Sampling problem (in the 
sense of sampling over an ensemble of large matrices us¬ 
ing as weight their permanents) is a particular case of a 
much more general kind of physical situations where the 
transition amplitudes between many-body Fock states 
built from two different single-particle basis sets are mea¬ 
sured. Within this general scenario, Boson Sampling re¬ 
quires the calculation of the many-body unitary oper¬ 
ator representing a linear, canonical transformation at 
the single-particle level. We have provided different ver¬ 
sions of the problem, obtained by expressing this tran¬ 
sition amplitudes in different intermediate basis like co¬ 
herent states and quadrature states of the field. Start¬ 
ing with these exact representations, we performed an 
asymptotic analysis valid in the limit of large occupa¬ 
tions and provide their semiclassical approximation in 


the spirit of coherent sums over solutions of a classical 
boundary problem. Along the way, we have derived an ex¬ 
act form of the many body transition amplitudes, equiv¬ 
alent to the calculation of permanents, and use it to de¬ 
rive exact results for the moments of the distribution of 
permanents over the Ginibre ensemble. Work on the ex¬ 
tension of our asymptotic analysis into the regime of low 
densities, where BS is expected to be hard, is currently 
under way. 

Acknowledgements. We thank M. Tichy, A. Buchleitner, 
J. Kuipers and V. S. Shchesnovich for valuable discussions 
and three anonymous referees for their help in improving the 
paper. 

Key words. Boson sampling, semiclassical approxima¬ 
tions, canonical transformations, many-body scattering, non¬ 
interacting quantum fields 


References 

[1] C. K. Hong, Z. Y. Ou, and L. Mandel Phys. Rev. Lett. 
59 , 2044(1987). 

[2] E. Farhi and S. Gutmann Phys. Rev. A 58, 915 (1998). 

[3] A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, 
P. I. Mosley, E. Andersson, I. lex, and C. Silberhorn 
Phys. Rev. Lett. 104 , 050502 (2010). 

[4] A. Schreiber, A. Gabris, P. P. Rohde, K. Laiho, M. Ste- 
fanak, V. Potocek, C. Hamilton, I. lex, and C. Silber¬ 
horn Science 336, 55 (2012). 

[5] A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, 

I. lex, and C. Silberhorn Phys. Rev. Lett. 106 , 180403 
( 2011 ). 

[6] S. Aaronson and A. Arkhipov, The Computational 
Complexity of Linear Optics, in: Proceedings of the 
Forty-third Annual ACM Symposium on Theory of 
Computing,, STOC ’ll (ACM, New York, 2011). 

[7] M. A. Broome, A. Fedrizzi, S. Rahimi-Keshari, 

I. Dove, S. Aaronson, T. C. Ralph, and A. G. White 
Science 339 , 794 (2013). 

[8] I. B. Spring, B. ]. Metcalf, P. C. Humphreys, W. S. 
Kolthammer, X. M. fin, M. Barbieri, A. Datta, 

N. Thomas-Peter, N. K. Langford, D. Kundys, I. C. 
Gates, B. I. Smith, P. G. R. Smith, and I. A. Walmsley 
Science 339 , 798 (2013). 

[9] A. Crespi, R. Osellame, R. Ramponi, D. I. Brod, E. F. 
Galvao, N. Spagnolo, C. Vitelli, E. Maiorino, P. Mat- 
aloni, and F. Sciarrino Nature Photon 7, 545 (2013). 

[10] M. Tillmann, B. Dakic, R. Heilmann, S. Nolte, A. Sza- 
meit, and R Walther Nature Photon 7, 540 (2013). 

[11] M. C. Tichy f. Phys. B: Atomic, Molecular and Opti¬ 
cal Physics 47 , 103001 (2014). 

[12] I. W. Negele and H. Orland, Quantum Many-particle 
Systems (Westview Press, 1998). 


10 


Copyright line will be provided by the publisher 
























Annalen der Physik, August 28, 2015 


[13] K. B. Wolf, Geometric optics on phase space 
(Springer, Berlin, 2004). 

[14] W. Vogel and D. G. Welsch, Quantum Optics, 3rd 
edition (Wiley-VCH, 2006). 

[15] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Gryn- 
berg, Photons and atoms : introduction to quantum 
electrodynamics (Wiley, New York, 1997). 

[16] W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. I. 
Gillen, S. Foiling, L. Pollet, and M. Greiner Science 
329, 547 (2010). 

[17] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and 
D. S. Fisher Phys. Rev. B 40 , 546 (1989). 

[18] M. Lewenstein, A. Sanpera Trigueros, and V. Ahufin- 
ger, Ultracold atoms in optical lattices : simulating 
quantum many-body systems, 1st edition (Oxford 
Univ. Press, Oxford, 2012). 

[19] S. Inouye, R. A. M., J. Stenger, H. J. Miesner, M. S. K. 
D., and W. Ketterle Nature 392,151 (1998). 

[20] P. Courteille, R. S. Freeland, D. J. Heinzen, F. A. van 
Abeelen, and B. J. Verhaar Phys. Rev. Fett. 81, 69 
(1998). 

[21] Y. V Fyodorov Int. Math. Res. Notices 2006 (2006). 

[22] V S. Shchesnovich Int. J. Quantum Inf. 11,1350045 
(2013). 

[23] W. H. Miller J. Chem. Phys. 53, 1949 (1970). 

[24] W. H. Miller Adv. Chem. Phys 25, 69 (2007). 

[25] J. D. Urbina, J. Kuipers, Q. Hummel, and K. Richter, 
arXiv:1409.1558. 


Copyright line will be provided by the publisher 


11 


